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Pure 5(7(3) glue theories exhibit a deconfining phase transition at a nonzero temperature, T c . 
Using lattice measurements of the pressure, we develop a simple matrix model to describe the 
transition region, when T > T c . This model, which involves three parameters, is used to compute 
the behavior of the 't Hooft loop. There is a Higgs phase in this region, where off diagonal color 
modes are heavy, and diagonal modes are light. Lattice measurements of the latter suggests that 
the transition region is narrow, extending only to about ~ 1.2 T c . This is in stark contrast to lattice 
measurements of the renormalized Polyakov loop, which indicates a much wider width. The possible 
implications for the differences in heavy ion collisions between RHIC and the LHC are discussed. 
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Heavy ion collisions at the Relativistic Heavy Ion Col- 
lider (RHIC) have demonstrated a rich variety of unex- 
pected behavior [1]. Notably, in peripheral collisions the 
elliptical flow can only be described by nearly ideal hy- 
drodynamics, with a very small ratio between the shear 
viscosity, 77, and the entropy density, s. The differences 
between collisions at RHIC, and those which will soon be 
observed soon at the Large Hadron Collider (LHC), will 
be especially interesting: does the nearly ideal hydrody- 
namic behavior, observed at RHIC, persist at the much 
higher energies of the LHC? 

One approach to deconfinement exploits the analogy 
to M = 4 supersymmetric gauge theories: using the 
AdS/CFT correspondence, such theories are computable 
analytically in the limit of infinite coupling, for an infinite 
number of colors [2] . By introducing a potential for the 
dilation field, the behavior of the entropy density near 
the deconfining phase transition, at a temperature T c , 
can be fit from measurements on the lattice [3-6] . While 
the entropy density, s, decreases strongly as T — > T+ 
because it is related to Hawking radiation, in AdS/CFT 
models the ratio 77/s remains completely independent of 
temperature. This suggests that like RHIC, that colli- 
sions at the LHC should also be described by nearly ideal 
hydrodynamics; see, also, Ref. [7]. 

In this work we consider a very different approach to 
the deconfining phase transition. It assumes that the 



coupling is moderate even down to the transition tem- 
perature, T c [8]. We use an elementary matrix model, 
involving three parameters, to parametrize the behav- 
ior of the deconfining phase transition. A version of 
this model with one parameter was first proposed by 
Meisinger, Miller, and Ogilvie [9]. Similar models arise 
for theories in which one (or more) spatial directions are 
of femtoscale size [10-13]. 

The parameters of the model are fixed from lattice 
measurements of the pressure [14-17]. It then predicts 
how the 't Hooft loop [18-22] changes with temperature 
near T c , which we compare to the results of lattice sim- 
ulations [23, 24]. Further, the model predicts that for a 
range of temperatures above T c , there is a Higgs phase, 
where correlation functions of electric fields are a mixture 
of heavy and light modes, from fields which are off diago- 
nal, and diagonal, in color, respectively. This may help to 
understand the results of lattice simulations [15, 25-27], 
which are otherwise somewhat puzzling. 

The most direct prediction of our model is for the ex- 
pectation value of the Polyakov loop. For the pure glue 
SU(3) theory, lattice simulations find that the (renormal- 
ized) Polyakov loop vanishes below T c , jumps to ~ 0.4 
at TJJ, and then rises with T, until it is approximately 
constant above ~ 4.0 T c [28-30]. This represents confine- 
ment below T c , a complete Quark Gluon Plasma (QGP) 
at high temperature, and a "semi" -QGP in between [31- 
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34] . Physically, there is no ionization of color in the con- 
fined phase, total ionization in the complete QGP, and 
only partial ionization in the semi-QGP [34]. (While we 
discuss a purely gluonic plasma, we adopt the common 
term QGP.) 

The principal thrust of this paper is that from indirect 
measurements on the lattice, we suggest that the width 
of the semi-QGP is much narrower than indicated by 
present results for the renormalized Polyakov loop: not 
up to ~ 4.0 T c , but only ~ 1.2 T c . We do not understand 
this discrepancy in detail, but suggest a possible reason 
later. This discrepancy is the reason why, having fit the 
parameters of our model from the pressure, we compute 
both the 't Hooft loop and gluon masses. 

While we treat the pure glue theory, our model can 
be extended to QCD, with dynamical quarks [17]. It is 
reasonable to assume that in QCD, the semi-QGP is like 
that of the pure glue theory, relatively narrow. We thus 
conclude by discussing the possible phenomenological im- 
plications of our results for heavy ion collisions. 

How confinement arises in our model can be under- 
stood by analogy. For a bosonic field, with energy E and 
chemical potential fi, the Bose-Einstein statistical distri- 
bution function is 
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Instead of taking \x to be real, as in ordinary thermody- 
namical systems, for the purposes of the analogy we take 
it to be purely imaginary, and define \x — lixTiq, where 
q is real. 

Doing so, n(E, q) is clearly a periodic function of q, 
invariant under q — > q + 1 . Thus we can choose to define 
q to lie within the range from — | to +i. 

Now assume that we integrate over q, with a distri- 
bution which is flat in q. Expand for large energy, so 
that the first term is the Boltzmann statistical distribu- 
tion function. Given the assumed distribution in g, the 
integral of this term vanishes, 
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Indeed, we can expand the Bose-Einstein distribution 
function term by term in powers of Boltzmann factors, 
e -E/T+27Ttq j-^gj. doing S0; the integral over each and ev- 
ery term obviously vanishes. The same is true for the 
Fermi-Dirac distribution function as well. 

Thus a flat distribution in q represents the confined 
phase. To represent a phase with partial deconfinement, 
one integrates over a limited region, say q : —qo — > +go , 
with go < \- Complete deconfinement occurs when one 
integrates over a distribution which is a delta-function in 

q- 

This example appears somewhat artificial. For a given 
g, the statistical distribution functions are complex val- 
ued, and so only integrals over q can possibly represent 



physical quantities. Indeed, the grand canonical ensem- 
ble is characterized by a fixed value for the chemical po- 
tential, and not by an integral over /x's. 

Nevertheless, precisely this mechanism arises for the 
deconfining phase transition in a SU(N) gauge theory. 
Consider the expansion about a background field for the 
time-like component of the vector potential, 



2nT 



q a S ab 
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a and b are colors indices, running from 1 . . . N. For 
nonzero g a 's, this background field acts like an imaginary 
chemical potential for the diagonal elements of the gauge 
group. Integration over the g a 's arises from imposing 
Gauss' law for those elements of the gauge group [19]. 

This background field generates a non-trivial expecta- 
tion for the Polyakov loop, £, which is the color trace of 
the thermal Wilson line, L: 



— trL 
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L = V exp ig 
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A dr 
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V represents path ordering, T is the temperature, and r 
the imaginary time, r : — > l/T. 

The lattice demonstrates that near T c , the expectation 
value of the Polyakov loop is not near one, and decreases 
as the temperature does. In such a region, the eigenvalues 
of the logarithm of the thermal Wilson line are nonzero. 
Taking an ansatz such as Eq. (3) is the simplest possible 
way to model this. We do not attempt to derive the 
distribution of these eigenvalues, but to guess that from 
lattice results. 

Since the gauge potential A is an element of SU(N), 
S^Li 9a = 0, modulo one, there are N — 1 independent 
g a 's. At infinite N, the g a 's form a continuum, and the 
example of Eq. (1) is exact; see, e.g., computations on 
a femtosphere at N — oo [13]. For two colors, we can 
choose the eigenvalues to be qi = —qi\ for three, q\ — 
-q 2 , and q 3 = 0. 

In the presence of the background field of Eq. (3), a 
potential for the g 's is generated at one loop order [18- 
22], 



V P t(q a ) = 



3 4^ 

a,b=l 



g^(l-ltol) - (N 2 -i)- 
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where q a b — qa — 96, defined modulo one. The minimum 
is at q a — 0, where — V p t(0) is the pressure for an ideal 
gas of gluons. 

The potential V p t{q a ) enters in computations of the 't 
Hooft loop. It is useful to consider deconfinement as a 
type of spin system. A pure SU(N) gauge theory has 
N degenerate vacua, where the thermal Wilson line L 
equals one of the TV roots of unity, 



(6) 
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j = . . . (N — 1). The usual vacuum, with j = and 
L = 1, corresponds to all q a = 0. A Z(N) vacua with 
j = 1 and L = e 2 "^! corresponds to N- 1 g a 's = 1/JV, 
and the remaining element = — 1 + 1/N. 

At high temperature in the complete QGP, the the- 
ory lies in one spin state, which we can choose to be 
j = 0. One can compute tunneling between two degen- 
erate vacua by constructing a box which is long in one 
spatial direction, with j = at one end, and j = 1 at 
the other. An interface between the two ordered states 
forms in the center of the box, with the interface tension 
between the two computable semi-classically, using the 
potential V p t(q a ) [18-22]. This interface is equivalent to 
a 't Hooft loop which wraps around the center of the box 
[20]. 

As the temperature decreases and T approaches T c , 
domains with j ^ form and grow in size. They become 
increasingly probable, until at T~ and below, as a spin 
system the vacuum is completely disordered, a sum over 
many spin domains. 

We want to add terms to the effective potential which 
model the transition to deconfinement. We could add 
perturbative corrections to V p t(q a ), which have been 
computed to ~ g 3 [21], but invariably they give q a = (or 
a Z(N) equivalent state) as the vacuum. With a com- 
plete theory, such as the monopole model of Liao and 
Shuryak [7] , this potential could be computed directly. 

We adopt a more modest approach, attempting to 
guess the form of the non-perturbative potential. We fit 
the coefficients which enter to lattice results for the pres- 
sure, and then use it to compute other quantities. The 
advantage of our approach is that we can compute quan- 
tities not just in, but near thermal equilibrium. Such 
quantities, like the shear viscosity [34], are much harder 
to extract on the lattice. 

Since the Polyakov loop is an order parameter for de- 
confinement, a natural guess is that the non-perturbative 
potential involves Z{N) invariant elements of the Lie 
group. The first such term is the adjoint loop [12, 31-36]. 
Instead, following the authors of Ref . [9] , and computa- 
tions of the 't Hooft loop [18-22], we write a potential 
which is a polynomial in the e^'s. 

There are several symmetries which any potential of 
the QaS must satisfy. It must be periodic in each q a , 
with q a q a + 1. It must also be invariant under Z(N) 
transformations, where N— 1 of the g a 's shift by 1/N, and 
the last element, by — 1 + 1/N. Lastly, if we interchange 
the ordering of the g a 's, we can change q a b —> qta = 
—q a b- These symmetries can be satisfied by constructing 
a potential as a function of q a b(l — qab)- 

We can still form an infinite number of terms by tying 
together the color indices in different ways; see, e.g., the 
examples at two [21] and three [13] loop order. We adopt 
the simplest approach, and take terms like those which 
arise at one loop order, Eq. (5), which involve a sum over 



one q ab : 

N 

V non {q a ) = T 2 T 2 (°i WC 1 - \lab\) 

a, 6=1 

+C2q 2 ab (l-\qab\) 2 + c 3 ) . (7) 

The model of Ref. [9] involves terms ~ c\ and C3; these 
are kept in fixed ratio, given by the second Bernoulli 
polynomial. Instead, we allow C\ and C3 to vary indepen- 
dently. This helps to avoid a pathology of the model of 
Ref. ([9]), where the pressure is negative at T c . 

We also introduce a term ~ C2q 2 b (l — q a b) 2 ', this is 
proportional to the perturbative term in Eq. (5), and is 
related to the fourth Bernoulli polynomial. 

We take all of the non-perturbative terms to be ~ T 2 , 
since lattice simulations indicate that in the pure glue 
theory, the leading corrections to terms ~ T 4 are ~ T 2 
[9, 10, 33, 37]. There is obviously no fundamental rea- 
son why other terms, such as those ~ 1, could not also 
appear. 

When the g^'s develop an expectation value, this rep- 
resents symmetry breaking for an adjoint scalar field, Aq, 
coupled to an SU(N) gauge field, the A^s [33]. As an 
adjoint scalar, though, there is no strict order parameter 
which distinguishes between the symmetric and broken 
phases. Thus there need not be a phase transition in go- 
ing from the symmetric phase, the complete QGP, to the 
"broken" phase, the semi-QGP. 

If there were such a phase transition, it would repre- 
sent a second transition, above T c , separate from that 
for deconfinement. While possible, in a pure SU(N) 
gauge theory lattice simulations only find evidence for 
one phase transition, at T c [14-17]. To avoid a phase 
transition between the complete and semi-QGP, it is es- 
sential that the non-perturbative potential has a term 
which is linear in the q a 's. Assume that the effective po- 
tential only involved terms such as ~ ?S,(1 — <Zaf>)™ for 
71 > 2. For small q a , these are of quadratic or higher 
order in the <7 Q 's, and of necessity, there would then be a 
phase transition when the q a 's developed a nonzero ex- 
pectation value. This transition might be of either first 
or second order, but there would be a phase transition. 
When ci ^ 0, though, a term linear in the q a 's ensures 
that there is no such phase transition. Instead, even for 
high temperature, there is always a small but non-zero 
expectation value for the g a 's, < q a > ~ 1/T 2 ; that is, 
the theory is always in a Higgsed phase. As we shall see 
however, this point is somewhat academic. For the pa- 
rameters relevant to two and three colors, the region in 
which Higgsing matters is very narrow. 

We remark that effective theories on the lattice often 
exhibit phases with broken symmetry [38] . The necessity 
of such a broken phase near T c does not seem to have 
been appreciated previously, though. 

To determine the parameters of the model we compare 
to lattice measurements of the pressure. For three colors, 
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this is illustrated in Fig. (1); for two colors, in Fig. (2). If 
p(T) is the pressure, and e(T) the energy density, then a 
more sensitive test of the fit is also to plot the interaction 
measure, A = e — 3p. Thus in each figure we plot p/T 4 , 
e/T 4 , and A/T 4 , both from the lattice, from Ref. ([14]) 
for two colors, and from Ref. ([15]) for three colors. 
The parameters of the fit are 

Cl =- .41488 ; c 2 = - 5.45957 ; c 3 = 0.21954 . (8) 

for three colors, and 



where 



ci = - 0.30267 ; c 2 = - 5.97440 ; c 3 = 0.18341 . 



(9) 



for two colors. 

While our model appears to involve three parameters, 
this is misleading. One parameter fixes the critical tem- 
perature, T c . A second is chosen so that the pressure 
vanishes at T c . Thus we really have only one free param- 
eter, which is tuned to fit the behavior of the pressure 
near T c . 

For two colors, our model exhibits unphysical behav- 
ior, as the energy density is negative below ~ 1% of T c . 
This might be corrected by adding further terms in the 
nonperturbative potential, such as higher Bernoulli poly- 
nomials. 

In any case, since we fix the pressure to vanish at T c , 
within our approximations the confined phase has van- 
ishing pressure. How to match to a more realistic de- 
scription of the confined phase is an important problem, 
which we defer for now. 

Given the effective Lagrangian, it is then straightfor- 
ward to compute the 't Hooft loop. In the complete QGP, 
the potential includes only the perturbative potential, 
V p t(q a ), Eq. (5); in the semi-QGP, it is a sum of this and 
the non-perturbative potential, V non (q a ). 

For two colors, as q2 = —qi there is only one indepen- 
dent direction, and it is direct to compute the tunneling 
path, and its associated action, analytically. The result 
for the 't Hooft loop is 



Z(g 2 ) = 1 - 0. 16459 g 2 (T) 



am = 4 " 2r2 £(a 2 ) (i - (r c /r) 2 ) 3/2 

{) 3v/6^(T) 1 - 0.908 (T c /T) 2 



(10) 



The factors involving T c /T are special to the semi-QGP, 
so that as T > T c , the result reduces to that in the 
complete QGP [18]. The function £(c/ 2 ) is the correction 
~ g 2 in the complete QGP; in plotting, we take g 2 (2irT) 
[8]. 

The 't Hooft loop vanishes at T c , as expected for a sec- 
ond order phase transition. From universality, the result 
in the Ising model is a(T) ~ (T - T c ) 2v , with 2v ~ 1.26; 
lattice results in a gauge theory [23] find 2v ~ 1.32 [23]. 
Our result, 2v — 1.5, is not too far off, as expected for a 
mean field theory. We note, however, that because of the 
term in the denominator, that the numerical agreement 
isn't close. This is presumably related to the unphys- 
ical behavior of the energy density near T c , mentioned 
previously. 

For three colors, in the semi-QGP the vacua is along A3 
(using the Gell-Mann notation), while the path for the 
't Hooft loop depends upon a change in As- The path 
was determined numerically, and lies along both A3 and 
\$. The action of the tunneling path was also determined 
numerically, and the result for the 't Hooft loop for three 
colors is illustrated in Fig. (3). (For N — 2, we take 
T c /A Ms = 1.31; for N = 3, 1.14. For the same value of 
T c /A A j s = 1.31, the results unexpectedly coincide.) 

Including £(g 2 ), the semi-classical computation of the 
't Hooft loop in the complete QGP agrees with lattice 
simulations above ~ 4.0 T c ; below that temperature, they 
agree with the result in the semi-QGP [24]. To obtain 
agreement, however, it is necessary to include the correc- 
tion £,(g 2 ); this is computed in the complete QGP, which 
is incorrect. Two things are required to compute £(<7 2 ) in 
the semi-QGP. First, the potential for constant q a needs 
to be computed to two loop order, expanding about the 
full potential, V p t(q a ) + V non (q a ). Second, corrections to 
one loop order need to be computed for the kinetic term. 
In the complete QGP this brings in new functions, the 
ip(q a ) [18]. Other functions could arise in the semi-QGP. 
For now, we defer these involved computations; since the 
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FIG. 1: Comparison of lattice results for SU(3) pure gauge 
to the model, for the pressure, energy density, and interaction 
measure. 
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FIG. 2: Comparison of lattice results for 5(7(2) pure gauge 
theory to the model, for the pressure, energy density, and 
interaction measure. 
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corrections ~ g 2 are large, ~ 50%, our results should be 
considered as tentative. 

Besides the 't Hooft loop, which is an interface ten- 
sion for an order-order interface at T > T c , the interface 
tension for the order-disorder interface, at T c , is also com- 
putable in our model. This only exists for a first order 
transition; for three colors, 



0.0258012 



T 2 



(11) 



It is necessary to compute the corrections ~ g 2 before 
comparing to lattice data, though. 

The parameters for three colors, Eq. (8), and two col- 
ors, Eq. (9), are similar; the difference is commensurate 
with a dependence on ~ 1/A^ 2 , with the coefficient of 
order one. We have then assumed that the parameters 
for three colors are close to those for higher N. We find 
reasonable agreement for the interaction measure to lat- 
tice results [16]. When N > 4, there is more than one 
't Hooft loop. Lattice simulations find that they obey 
Casimir scaling to good approximation [24]. We have not 
checked this explicitly, but suspect that in our model, 't 
Hooft loops respect Casimir scaling. 

The most novel prediction of our model is that there 
is a Higgs effect in the semi-QGP. This was noted first 
in Ref. [33], and in theories at a femtoscale [12]. To 
understand it, consider the quantum fluctuations about 
the background field of Eq. (3): 



((AD ab (x)(A^) ba (0)) 



d 3 p 
(2tt) 3 



n— — oo 



where A o is the quantum propagator 



(12) 



A 00 = 



(p) 2 +P 2 + ^W) ;po = 2nT{n + qa - qb) 



(13) 
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FIG. 3: The 't Hooft loop for SU(3) pure gauge theory: lattice 
data from Ref. [24], and GKA, Giovannangeli and Korthals 
Altes, Ref. [21], the semi-classical computation in the com- 
plete QGP, including corrections of ~ g 2 . In our model we 
show results for two and three colors, assuming that the cor- 
rections of ~ g 2 are identical in the complete and semi-QGP; 
see text. 



The shift in the energy, p = 2nTn —t 2irT(n + q a — qi,), is 
because we are expanding about a background field. The 
background field Aq acts upon quantum fluctuations like 
an adjoint Higgs field. Because the field is diagonal in 
color, Eq. (3), diagonal fields do not feel the background 
field. Thus for diagonal fields, the only mass they develop 
is the Debye mass, tod- This is of order ~ gT times a 
function of the q a 's [34]. In contrast, off diagonal fields 
have non-trivial commutators with a diagonal field, and 
so they develop "masses" which are large, ~ 2wT(n + 
Qa - Qb)- 

We illustrate this in Fig. (4) for three colors. The 
masses of the two diagonal gluons are equal, and decrease 
as T — > T+. There arc two types of off-diagonal gluons: 
four with \a — b\ — 1, and two with \a — b\ = 2. The 
splitting of the masses is evident only close to T c , for 
T < 1.2 T c . 

We do not plot lattice data, because it is somewhat 
contradictory. Lattice measurements of a gauge invariant 
quantity, the two point function between Polyakov loops, 
shows that the associated mass decreases as T — > T+ [15]. 
In contrast, the two point function of gluons indicate that 
the gauge dependent mass increases as T — > T+ [27]. 
Clearly it would be best to reanalyze the lattice data 
with a Higgsed propagator in the effective theory, Eq. 
(13), with its characteristic combination of modes whose 
masses both increase and decrease. 

The static, spatial gluon fields, the Aj, also undergo a 
Higgs effect. This happens as well in a monopole gas [7]. 

We have not discussed the most obvious application of 
our model: the computation of the Polyakov loop. We 
plot this quantity, and the lattice results, for three colors 
in Fig. (5). A direct comparison of the two is somewhat 
misleading. We have not computed perturbative correc- 
tions to the Polyakov loop, which enter at ~ g 3 [39]. 
This contribution is positive, and will increase the re- 
sult. Nevertheless, while the two coincide at T c — which 
is presumably coincidence — they immediately diverge 
from one another. From Fig. (5), in our model the loop 
quickly goes up to a constant value by ~ 1.2 T c ; this is 
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FIG. 4: Gluon masses, m a t, / (gT) , for 5(7(3) pure gauge the- 
ory: diagonal gluons, with a = 6, are light; while there are 
two off diagonal gluons with heavy masses, for \a — b\ = 1 and 
la -61 = 2. 
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very different from lattice measurements, for which it is 
not constant until a much higher temperature, ~ 4.0 T c 
[28-30]. 

If our model is correct, why does the value of the 
Polyakov loop, computed from our model, differ so signif- 
icantly from lattice measurements of the (renormalized) 
Polyakov loop? There is an ambiguity associated with 
the renormalized Polyakov loop, from the zero point en- 
ergy. In Ref. [34] we argued that perturbatively, the zero 
point energy vanishes for a straight Polyakov loop. This 
argument fails for a "smeared" loop (see, e.g., the ap- 
pendix of Ref. [29]). If so, then the effects of smearing 
must be very dramatic. 

We comment that a similar rapid growth in the 
Polyakov loop is found in solutions of the Schwinger- 
Dyson equations [40]. Our results do not coincide nu- 
merically, though. 

To understand our results better, consider first the 
limit of intermediate temperature: say above ~ 1.5 T c , 
and up to ~ 4.0 T c . In this limit the <7 a 's are small. 
As noted before, for small q the dominant terms are 
~ c\T 2 q in the nonperturbative potential, Eq. (7), and 
that ~ T A q 2 in the perturbative potential, Eq. (5). Bal- 
ancing these two gives (q a ) ~ —ci/n 2 (T c /T) 2 ; for three 
colors, (q a ) ~ 0.04 (T c /T) 2 . Thus while the theory is 
nominally always in a "Higgsed" phase, as a practical 
matter this effect is numerically miniscule. Further, since 
the loop involves the cosine of the g a 's, asymptotically 
the deviation of the loop from unity is even smaller: 
1 - £~ q 2 ~ (T c /T) 4 . 

At intermediate temperatures, what is the origin of 
the term ~ T 2 in the pressure; or equivalently, the be- 
havior of the interaction measure, A/T 4 ~ 1/T 2 ? If the 
q ~ 1/T 2 , for small q terms ~ T 4 q 2 and ~ T 2 q are of 
order one. Then the only contribution to a term ~ T 2 
in the pressure is from the g-iraiependent terms in the 
nonperturbative potential. In our model, this is a single 
term, ~ c 3 T 2 T 2 . 

On the other hand, for intermediate temperatures, 
even if the pressure does not probe the other terms in 
the potential, the 't Hooft loop does. By measuring a 
't Hooft loop, the altered boundary conditions force the 
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FIG. 5: The Polyakov loop for a SU(S) pure gauge theory 
from lattice simulations [30] and in our model. 



system to probe the g-dependence of the potential in a 
nontrivial way. In the present model, these are deter- 
mined by the coefficients for c\ and c 2 . (As a constant 
term, the value of C3 doesn't matter.) In this way the 
lattice measurements of the 't Hooft loop [23, 24] are an 
absolutely essential constraint on our model. 

In contrast, near T c all of the parameters of the model 
matter, contribute to both the behavior of the pressure 
and the value of the <? a 's. By comparison, consider the 
model of Ref. ([9]), where c 2 = 0. We take two colors, 
since then the expressions are algebraically simple. From 
Ref. ([9]), with q = q x = -q 2 , 

«~- i( 1 -i/ 1 -(t)') ■ < 14) 

When c 2 = 0, though, the interaction measure is much 
broader than indicated by the lattice data. 
In the present model, with c\,c 2 7^ 0, 

?= i( 1_ l/ 1 + ^(^)(w) ■ (15) 

The values of c± and c 2 are dictated by fitting the pres- 
sure, or more sensitively, the interaction measure. By 
making c 2 increasingly negative from c 2 = 0, one finds 
that the peak in the interaction measure becomes increas- 
ingly sharp. Since the loop vanishes at T c for two colors, 
ci is then adjusted so that q = 1/4 at T = T c . 

Tuning c 2 and c\ in this way, one finds that as the 
peak in the interaction measure sharpens, that the re- 
gion in which q is nonzero narrows as well. This is due 
to the particular values of our fit: near T c , the term 
~ c 2 T 2 T 2 is not only large, but approximately cancels 
the similar term, with coefficient (27r 2 /3) T 4 , in the per- 
turbative potential. Requiring q = 1/4 at T = T c fixes 
c 2 + 27r 2 /3 = — 2ci, so that |ci| is small, |ci| <C \c 2 \. From 
Eq. (15), the combination ci/(c 2 + (2tt 2 /3)(T c /T) 2 ) en- 
ters into q, and implies that it is much sharper than the 
corresponding factor, (T/T c ) 2 for c 2 = 0, Eq. (14). A 
similar cancellation happens for three colors. 

We suggest that this reflects real physics: the sharp- 
ness of the interaction measure reflects the narrowness of 
the region in which the loop deviates significantly from 
unity. For the model with c 2 = 0, Eq. (14), q al- 
ways exceeds the corresponding value in our model, with 
C2 / 0. This is also seen for T ^> T c : when c 2 = 0, 
q ~ 0.125 (T c /T) 2 ; with c 2 ^ 0, q ~ 0.011 (T c /T) 2 . 

It is also worth contrasting our results with those in 
a Polyakov loop model. Consider a theory which only 
involves the Polyakov loop of Eq. (4), 

V ef f(l)= (-|l^| 2 + J(K| 2 ) 2 )& 4 T 4 ; (16) 

see, e.g., Eq. (2) of Ref. [32], and Polyakov Nambu- 
Jona-Lasino (PNJL) models [35]. For three colors, the 
Z(3) symmetry also allows a cubic term, ~ £ 3 + (£*) 3 , 
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but its addition would only complicate the algebra, and 
not our qualitative conclusion. There is no cubic term 
for two colors. 

The minimum of the potential is £q = \fb-2, which we 
choose to be real. As it is related to the pressure of 
an ideal gas of gluons, we assume that the coefficient 64 
is independent of the temperature, and that only 62, or 
equivalently to, depends upon T. The pressure and the 
interaction measure are then 

-^ = it; ^ = Ht 9 ^. (17) 

Pideal ^Pideal OT 

Consider first intermediate temperature, where the ex- 
pectation value of the loop is near one. To obtain a 
term ~ T 2 in the pressure, in a loop model it is nec- 
essary to assume that the loop deviates from one as 
1 — £0 ~ (T c /T) 2 + . . .. This is contrast to our matrix 
model, where the deviation of the loop is ~ l/T" 4 , but 
there is still a term ~ T 2 in the pressure. 

Near T c , from Eq. (17) a peak in the interaction mea- 
sure corresponds to a rapid change in £q. This is similar 
to what we find in a matrix model. For I to decrease as 
T — > T+, so must in Eq. (16). This is proportional 
to the mass of the £ field, and so the t mass decreases, like 
that of the diagonal modes in the matrix model, Fig. (4) ; 
it is in contrast to the masses of the off-diagonal modes, 
which increase. 

In Ref. [32] and other loop models [35], in order to fit 
the pressure the temperature dependence of the pressure, 
62 must have a complicated form. In our matrix model, 
the coefficients are just - T 4 , Eq. (5), and - T 2 , Eq. 
(7). In a mean field theory such as this, simplicity is a 
virtue. 

Lastly, the splitting of gluon masses near T c is special 
to a matrix model, as a Higgs effect for the adjoint scalar 
Aq field. See, e.g., the loop model of Ref. [10], where the 
splitting of masses does not occur. 

Our analysis is a preliminary first step. In deriv- 
ing our results, we balance the perturbative potential, 
Vpi(g , ), against the non-perturbative potential, V non (q a ). 
In powers of g 2 , the perturbative potential is of order one, 
so implicitly we have assumed that the non-perturbative 
is as well. Since the non-perturbative potential repre- 
sents a resummation of effects to all orders, this is a 
strong assumption. Nevertheless, it allows us to envis- 
age computing to higher order in g 2 . Corrections at least 
to order ~ g 2 and ~ g 3 are needed in order to make 
a serious comparison to lattice data. This also requires 
precise lattice data, close to the continuum limit, not just 
for the pressure, but also for the 't Hooft loop and gluon 
masses. 

There are several formal questions raised by our analy- 
sis. The parameters of effective theories can be computed 
from lattice simulations [36] ; doing so for elements of the 
Lie algebra, instead of for elements of the Lie group, may 
be useful. It is also necessary to extend the analysis of 
Hard Thermal Loops in the complete QGP to the semi- 
QGP. This is equivalent to understanding the analytic 



continuation of the thermal Wilson line from imaginary 
to real time. 

To compare with QCD it is necessary to include the ef- 
fects of dynamical quarks. It will be especially interesting 
to see if, upon adding the effects of quarks to the pertur- 
bative potential, V p t(q a ), whether the thermodynamics 
[17], and the Debye mass, are reproduced using the same 
parameters for the non-perturbative potential, V non (q a ), 
in the pure glue theory. (With dynamical quarks, the 't 
Hooft loop does not exist as an order parameter.) 

Without detailed computation, we assume that a nar- 
row width for the semi-QGP in the pure glue theory im- 
plies the same for QCD. We thus conclude with some 
speculations for the phenomenology of heavy ion colli- 
sions. 

If RHIC probes to some temperature in the QGP, then 
LHC may probe to a temperature approximately ~ 50% 
higher. If the AdS/CFT correspondence holds for QCD, 
then results at the LHC must mimic those at RHIC. With 
the present analysis, the picture is rather more compli- 
cated. 

We assume, for the sake of argument, that RHIC 
probes only to a temperature in the semi-QGP, very near 
T c . Then the LHC begins at a temperature well in the 
complete QGP. Any conclusions are tempered by the fact 
that even if the LHC starts at a higher temperature, as 
it cools it must traverse through the semi-QGP. 

In the semi-QGP, the ratio of -q/s decreases as the 
square of the Polyakov loop as T — >• T c + ; this is true both 
in the pure glue theory, and with dynamical quarks [34] . 
Conversely, then, rj/s increases as the temperature goes 
up from T c . This is in sharp contrast to models based 
upon the AdS/CFT correspondence, where rj/s is con- 
stant [3-5]. Unfortunately, a computation beyond lead- 
ing logarithmic order is required to compute the precise 
dependence of rj/s with temperature. 

If the shear viscosity increases strongly from T c , and 
the system is in thermal equilibrium, then an increased 
shear viscosity should lead to an increase in particle mul- 
tiplicity, and a decrease in the elliptical flow, over the 
results expected from a (nearly) ideal gas. If the shear 
viscosity increases significantly, though, a hydrodynamic 
description could easily break down. 

It is also possible that the temperature dependence 
of 77 /s is weak; if so, then the particle multiplicity and 
elliptical flow at LHC should be similar to that expected 
by an extrapolation from the results at RHIC. There are 
then other ways to probe the effects of the semi-QGP. 

Consider, for example, energy loss, which is controlled 
by a parameter q. In the complete QGP, q ~ T 3 , or 
equivalently, the entropy density, s. In kinetic theory, 
q/s and s/rj are each proportional to a cross section, 
so one expects that a minimum in rj/s corresponds to 
a maximum in the energy loss, q/s [7, 41]. Following the 
methods of Ref. [34], the energy loss of a quark can be 
computed in the semi-QGP; as T — > T c + , it vanishes lin- 
early in the Polyakov loop. Thus in the semi-QGP, both 
T]/s and q decrease asT-> T+; the difference from Refs. 
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[7, 41] is because the kinetic theory for the semi-QGP is 
in the presence of a background Aq field. As the tem- 
perature increases from T c , then, excluding the obvious 
dependence upon the entropy, the energy loss is larger 
in the complete QGP than in the semi-QGP. As with 
the shear viscosity, determining the precise dependence 
upon temperature requires computation beyond leading 
logarithmic order. 

There is also a qualitatively new phenomenon in the 
semi-QGP: besides energy loss, the propagation of a col- 
ored field is suppressed by the background A field [34]. 
This suppression is universal, independent of the mass or 
momentum of the colored field. A complete analysis need 
incorporate this universal suppression as well as energy 
loss. 

Lastly, we note that given the modified propagator of 
the semi-QGP, Eq. (13), there are also significant modi- 
fications to the heavy quark potential [42] . This can also 
be compared to lattice data, which we defer for now. 

In the end, our speculations will soon be rendered moot 
by the wealth of results which will flow from heavy ion 
collisions at the LHC. The present approach is based 
upon constructing an effective theory from the results 
of lattice simulations; not just of the pressure, but quan- 



tities such as the 't Hooft loop and screening masses. 
This can then be tested against predictions from the 
AdS/CFT correspondence [2-6] and other models [7]. 
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